Attosecond probing of instantaneous AC Stark shifts in helium atoms 
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Based on numerical solutions of the time-dependent Schrodinger equation for either one or two 
active electrons, we propose a method for observing instantaneous level shifts in an oscillating strong 
infrared (IR) field in time, using a single tunable attosecond pulse to probe excited states of the 
perturbed atom. The ionization probability in the combined fields depends on both, the frequency 
of the attosecond pulse and the time delay between both pulses, since the IR field shifts excited 
energy levels into and out of resonance with the attosecond probe pulse. We show that this method 
(i) allows the detection of instantaneous atomic energy gaps with sub-laser-cycle time resolution 
and (ii) can be applied as an ultrafast gate for more complex processes such as non-sequential 
double- ionization. 

PACS numbers: 42.50.Hz, 32.80.Rm 



The energetic shift of atomic levels in external elec- 
tric fields is a well-known phenomenon and usually re- 
ferred to as "Stark shift" . For static fields that are much 
weaker than intra-atomic Coulomb fields, Stark shifts can 
be calculated using perturbation theory [1]. For oscil- 
lating external fields in the optical and near-IR range, 
perturbation theory breaks down at intensities of 10^^ 
W/cm^ [2], orders of magnitudes below the peak inten- 
sities available in state-of-the-art ultrashort laser labo- 
ratories. If such strong external fields are maintained 
over many optical cycles, cycle- averaged level shifts can 
be evaluated, e.g. by exploiting the quasi periodicity of 
the external field using the non-perturbative Floquet the- 
ory [3-5]. However, for the recently developed strong few- 
cycle IR laser pulses [il-fiot. atomic level shifts are non- 
perturbative in nature and also render the continuum- 
wave Floquet picture inappHcable. 

Modern pump-probe experiments combine extended 
ultraviolet (XUV) attosecond pulses of sub-IR-cycle pulse 
lengths (1 as = 10~^^ s) with phase-coherent IR laser 
pulses to observe electronic dynamics in atoms, molecules 
and solids [lll-[l6|. The role of laser-dressed highly ex- 
cited energy levels in atomic excitation and ionization 
has been studied recently using attosecond technology 
17 , [ist . In this Letter, we show that the attosecond 
pump-probe technique should also enable the measure- 
ment of instantaneous level shifts of low lying bound 
atomic states in alternating optical electric fields. Fur- 
ther, we demonstrate how the control of instantaneous 
level shifts can be exploited to gate other strong-field 
phenomena, such as non-sequential double ionization 



(NSDI). 

To this end, we simulate an XUV pump - IR probe 
scenario. We choose IR laser fields with negligible distor- 
tion of the ground state of He, that are, however, strong 
enough to couple low excited and continuous states, in- 
ducing noticeable level splitting, shift, and decay. For 
the XUV pulses we fix the number of cycles and vary the 
central frequency of the pulse. Key to our investigation is 
the observation that, for a given central frequency uosa of 
the single attosecond (SA) pulse and depending on the 
delay At between pump and probe pulse, the IR pulse 
may shift low-lying bound states into or out of resonance 
with one-photon excitations from the He ground state. 
The excited atom may then be easily ionized by the IR 
pulse. If the SA pulse is applied while the instantaneous 
level energies are off (in) resonant with uosa-, less (more) 
excitation and thus less (more) ionization out of excited 
states is expected to occur. This suggests that detection 
of the ionization probability as a function of uosa and At 
can be used to track the instantaneous Stark shifts. 

We will analyze this strategy by modeling He in the so- 
called single-active-electron approximation (SAEA), in 
which the time-dependent Schrodinger equation (TDSE) 
in velocity gauge reads (unless indicated otherwise, we 
use Hartree atomic units, e = m = h = 1): 
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V{z,p) ^{z,p;t), 
(1) 

where z and p are cylindrical electronic coordinates par- 
allel or perpendicular to the laser-polarization axis, re- 
spectively, and Pz and Pp are the corresponding conju- 
gate momentum operators. A{t) = —c J dtE{t) is the 
vector potential associated with the XUV and IR laser 
fields and c the speed of light. V{z^p) models electronic 
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FIG. 1: (Color online) Energy levels for He, calculated in 
SAEA without (black line) and with static external electric 
fields with intensities of 10^^ W/cm^ (red dash-dotted line) 
and lO^^W/cm^ (blue dashed line). 



correlation in terms of screening of the nuclear Coulomb 
potential by the passive electron [l9| . 

For later reference, we first compute energy- level 
shifts in the static field Est in the length gauge. 
Fig. [T] shows the spectrum S{E) of He atoms with- 
out external field and in static fields with intensities 
of 10^^ and lO^'^W/cm^, calculated by Fourier trans- 
formation of the autocorrelation function, S{E) = 
I /dt(V^(t)|V^(0))e-^^^P [20|. At these intensities the 
ground state is nearly unaffected by the external static 
field, while the excited 2s and 2p states are shifted to 
lower and higher energies, respectively [2l|. The level 
spacing e[^^2^^'^^ between the field-shifted Is and 2s levels 

decreases with increasing Egt^ while £^^^*2p^^ increases. 

Exposed to the combined electric field of a SA pulse 
and a delayed IR pulse Eiji{t) 



E{t) = sin(cj5A^) exp 



-^EiR ^m{ujiRt) cos^ 
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with EjR = if \t\ > Tir/ 2^ level shifts induced 
by EiR{t) can be probed on an attosecond time scale. 
We choose Gaussian SA pulses with an intensity of 
2 X 10^^ W/cm^ a pulse duration tsa (FWHM) of two 
XUV cycles, and variable uosa^ and a cos^-shaped IR laser 
pulse with a central wavelength of 800 nm, an intensity 
of 3 X 10^^ W/cm^, and a pulse duration tjr of four IR 
cycles. We neglect the spatial intensity profile of the IR 
laser pulse since SA pulses can be made to only overlap 
with the spatial center of the IR pulse [lH, |23| . 

We calculate the ionization probability of He in the 
oscillating field, Eq. ([2]), by propagating Eq. ([T]) on a nu- 
merical grid with equidistant spacings = Ap = 0.3 
and time steps 5t = 0.05. The spatial grid includes 2000 




FIG. 2: (Color online) Ionization probabilities (logarithmic 
color/grey scale) of He calculated in SAEA as a function of 
the center frequency ujsa of the SA pulse and time delay At 
between the SA and IR laser pulses in units of the IR laser 
period Tir. Superimposed dashed and dotted curves show the 
quasi-static level spacings ^^5^2^^^* ^^s^2p^^ respectively. 
Laser parameters are given in the text. 



(400) points and covers the range from —300 to 300 (0 
to 120) along the z {p) axis. We use a cos^/^ mask- 
ing function to suppress reflections from the grid bound- 
aries [2^ and determine the single-ionization probability 
from the accumulated outgoing electronic current into 
the absorber region. The wave function is propagated 
until the ionization probability becomes stabilized. 

The single ionization probability Ps is shown in Fig. O 
For At < the SA pulse precedes the center of the IR 
laser pulse. Accounting for the large bandwidth of the 
SA pulse, enhanced ionization may occur, even for usa 
below the ionization threshold, via excitation to a IR 
laser-dressed bound state followed by ionization in the 
remaining IR field. Single ionization is thus mediated 
by either the high-energy side of the XUV-pulse spec- 
trum that extends above the ionization continuum or in 
a two-step process via transiently excited atoms. Since 
Psi^SA^ At) behaves differently for different ujsa^ we an- 
alyze separately three usa intervals. In order to support 
our interpretations we present in Fig. [3] (normalized) ion- 
ization probabilities as functions of (a) At and (b) ujsa 
fixing the other parameter to different values. 

(1) For ujsA < 0.6 and negligible pulse-overlap (|At| > 
l.bTjR^ where Tir is the period of the IR field), there 
is no enhancement of the ionization probability. How- 
ever, Ps is enhanced and oscillates in the pulse-overlap 
region (|At| < 1.5Tir) as a function of At with period 
Tir/ 2. Ps reaches maxima when the SA pulse is applied 
at maximal IR electric-field magnitudes \Eir\ (Fig. [31(a), 
red-dashed line). We explain this behavior as follows. 
For I At I > 1.5T/i^, the energy gap between the (essen- 
tially) field- free Is and 2p level, E^l^2p — 0-75, is larger 
than UJSA^ and the spectral width of the SA pulse does 
not overlap with the 2p level. Thus excitation by the SA 
pulse is negligible. In contrast, for |At| < l.bTiR the 
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energy gaps £^^^^2!^^^ (^) ^u,2p^^ (^) static exter- 

nal field Est = Eiji{t), where EiR{t) is the instantaneous 
electric field of the IR-laser pulse at time t (cf. Fig. [1]). 
The results are drawn as dashed and dotted curves in 
Fig. [2] and match the contours of the ionization probabil- 
ity. 

For a quantitative investigation, we present the ion- 
ization probability Ps in Fig. [3](b) as a function of ujsa 
at three different delays. At At = —Tir/A (correspond- 
ing to a maximum of \EiR{t)\^ red dashed line) the in- 
stantaneous energy spacings can be identified by the two 
local extrema in Ps. The positions of these maxima as 
a function of ujsa match the corresponding quasi-static 
energy gaps ^^^^(t) and Ef^^^;'\t) for Est = \EiR{t)\ 
at t = -Tir/A. 

In contrast, application of the SA pulse at a zero of 
the IR field (At = 0, black solid line) leads to just one 
maximum in the ionization probability. The position of 
this maximum as a function of ujsa is slightly shifted 
to lower frequencies as compared to the field-free energy 
spacing £^^5^2^ (blue dash-dotted line for At = —2Tir in 
Fig. [3Kb)). This is due to the fact that the SA pulse has 
a finite non-zero temporal width and probes the (field- 
dressed) 2s and 2p levels near a zero of the IR field as 
well. 

Understanding of the ultrafast shift of the electronic 
levels can be useful to control electronic transitions with 
a high temporal resolution. Single ionization of an atom 
(or a molecule) is the initial step in many strong-field 
processes, such as higher-order harmonic and attosecond- 
pulse generation, NSDI, or molecular dissociation. Our 
findings above suggest that these processes can be gated 
in a controlled way on an attosecond time scale. We 
demonstrate the potential of this gating technique for 
NSDI of He by solving the TDSE 
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FIG. 3: (Color online) (a) Ionization probabilities, normalized 
to their respective maxima, as a function of At for three values 
of uJsA- 0.5 a.u. (red dashed curve), 0.75 a.u. (blue dash- 
dotted curve), and 1.1 a.u. (black solid curve), (b) Ionization 
probabilities as a function of ujsa for three values of At/Tm: 
—2 (blue dash-dotted line), —0.25 (red dashed line), and 
(black solid line). 



dressing by the IR-laser field repeatedly shifts 
into resonance with larger values of ojsa^ allowing for ef- 
ficient excitation followed by ionization in the IR field. 
The Tir/2 oscillation of Ps indicates that £^^^^23^^^^ 
pends on the instantaneous IR intensity. 

(2) If ujsA approximately matches the field-free energy 
gap E^^l^p and At < —1.5Tir, excitation to the 2p state 
by the SA pulse followed by subsequent ionization via the 
IR pulse is very likely. In contrast, the IR pulse by itself is 
not able to excite and leaves the atom in its ground state. 
Thus, for At > 1.5T/i^, the subsequent not-IR-assisted 
SA pulse can only excite, but not ionize, the atom. In 
the overlap region (|At| < 1.5T/i^) a small- amplitude IR- 
laser-induced oscillation of Ps occurs. Maxima in Ps 
occur at vanishing instantaneous IR intensity (Fig. [3l[a), 
blue dashed-dotted line), since the energy gap E[^^^2p^^^^ 
increases with field strength which leads to a mismatch 
with USA' In Fig. [3] the blue dashed-dotted line and 
red dashed line have opposite carrier phases, since the 
corresponding energy gaps are shifted into and out of 
resonance with the SA frequency, asynchronously. 

(3) For ujsA > 0.8 the spectrally broad SA pulse can di- 
rectly ionize He from the ground state, leading to notice- 
able ionization probabilities even without pulse overlap 
(Fig. [3l^a), black solid line). The larger ionization prob- 
ability for At < — 1.5T/i^ is due to dominant SA-pulse 
excitation followed by IR ionization, while this two-step 
mechanism does not apply for At > 1.5T/i^. In the over- 
lap region |At| < 1.5T/i^, the red dashed and black solid 
lines are in phase due to synchronous shifts into and out 

r -ji 7-,(dressed) 7-,(dressed) ,. -, 

01 resonance with i^i5^2s ^is,2p 5 respectively. 

Further evidence for the observation of Ps{ujsa^ ^t) 
allowing the detection of instantaneous level shifts is ob- 
tained by comparison with the instantaneous quasi-static 
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within a correlated two-electron model [2f 
field-free Hamitonian 
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where 
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{Z H- (-l)%/2)2 + p74. In this model 



the electronic center-of-mass coordinate Z is constrained 
along the polarization direction, z and p represent the 
coordinates of the (unconstrained) relative motion of the 
electrons parallel and perpendicular to the polarization 
axis, respectively. P^, Pz and Pp are the corresponding 
conjugate momentum operators. The soft core parameter 
s = 0.135. We solve Eq. (j3]) on a three-dimensional nu- 
merical grid with equidistant spacing AZ = Az = Ap = 
0.3 and 600, 1200, and 200 grid points along the Z, 
and p axis, respectively. To quantify the probabilities 
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FIG. 4: Probabilities for single ionization (SI) (a) and double 
ionization (DI) (b) as a function of the time delay between 
SA and IR pulses. The SA pulse has a central frequency of 
(^SA = 0.76 a.u. and a peak intensity of 2 x lO^^W/cm^. 
The IR pulse has a central wavelength of 800 nm and a peak 
intensity of 3 x lO^^W/cm^. The simulation results (circles) 
are interpolated by Unes. 

of single and double ionization, we partition the grid as 
ri < 12 and r2 < 12 (neutral helium), < 6 and rj > 12 
7^ J = 1? 2, He+) and the complementary space (He^+) 

(c.f. H). 

Fig. m^a) and (b) show the single and double ionization 
probabilities as functions of At between the SA and IR- 
laser electric fields ([2]). The SA (IR) pulse has a central 
wavelength of 60 (800) nm and peak intensity 2 x 10^^ 
(3x 10^^) W/cm^. The value of us a = 0.76 is chosen such 
that it is repeatedly in resonance with the field-dressed 
1^25 level. Note that e[^^^2] = 1 in this two-electron 
model. Therefore, FigHl (a) closely resembles in shape 
the USA < 0.6 result in the SAEA calculation in Fig. 
[3l^a), and the single- ionization-probability maxima occur 
when the SA pulse coincides with the maxima of \EiR{t)\. 

Once released, an electron can be accelerated and 
driven back to the core by the IR laser field, causing 
non-sequential rescattering ionization of the ionic core 
after about 2/3 Tir [26|. Hence, the single and double- 
ionization probabilities oscillate with the same period, 
Tir/2 (Fig. [4]), with maxima at almost identical de- 
lays. The double ionization probability almost vanishes 
for At > 0, since the decreasing IR laser intensity does 
not transfer sufficient energy to the rescattering electron 
to ionize the He+ core. If the SA appears much earlier 
than the IR pulse (At < — 1.5T/i^), the SA pulse first 
excites the helium atom, which is later singly and doubly 
ionized, once nearly maximal IR intensities are reached 

In conclusion, our simulations indicate that by apply- 



ing a SA pulse instantaneous energetic shifts in singly 
excited He induced by a strong few-cycle IR laser pulse 
are mapped onto oscillations in delay-depended single 
ionization probabilities. We interpret these oscillations 
in terms of a two-step process, where excitation in the 
SA pulse is followed by efficient ionization out of ex- 
cited states in the IR field. The excited-state population 
depends on the instantaneous energy gaps between the 
ground and excited states in the IR-laser field. Knowl- 
edge of the instantaneous energy levels in the strong IR 
field may lead to new schemes for the coherent control of 
NSDI, high harmonic generation, and molecular dissoci- 
ation. 
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